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Abstract 

Non-linear electron-screening corrections of stellar nuclear fusion rates are 
calculated analytically in the framework of the Debye-Hiickel model and com- 
pared with the respective ones of Salpeter's weak screening approximation. 
In typical solar conditions, the deviation from Salpeter's screening factor is 
less than one percent, while for hotter stars such corrections turn out to be 
of the order of one percent only over the limits of the Debye-Hiickel model. 
Moreover, an investigation of the impact of such non-linear screening effects 
on the solar neutrino fluxes yields insignificant corrections for both the pp and 
CNO chain reactions. 
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I. INTRODUCTION 



Stellar nuclear reactions rates are influenced by the electron cloud that screens fusing 
nuclei from each other. In recent years, there has been a wide interest of the actual effect of 
such screening corrections on reaction rates. That interest stems mainly from the fact that 
any variation of the fusion reaction rates reflects on the solar neutrino fluxes. Therefore the 
investigation of the screening effect has produced a number of papers, most of which come 
to the conclusion that the solution of the neutrino problem cannot be found in the screening 
effect, though it has been argued that the discrepancy between theory and experiment might 
be reduced by means of a suitable screening correction fl[]. Investigations pointing out the 
inability of screening corrections to reconcile theory and experiment are for example that of 
Ricci et all [[| and that of Gruzinov and Bahcall 0|. The former modifying Mitler's model || 
showed that possible uncertainties due to screening are too small to solve the solar neutrino 
problem while the latter concluded that Salpeter's weak screening formula is adequate for 
most solar fusion reactions, thus proving that screening can only modify solar reaction rates 
and neutrino fluxes by a few percent. 

However, as there is an exhaustive effort for higher precision in the theoretical as well as 
experimental calculation of neutrino fluxes, even corrections of a few percent have become 
significant || 

Most of the studies in that field depart from Salpeter's weak screening formalism [H] , and 
consider various corrections such as those arising from vacuum polarization [0] , electron den- 
sity distribution || etc. However, even authors who attempt to go beyond the linear regime 
make assumptions that lead to inaccuracies. For instance, a popular oversimplification is 
to consider that the electron density around the nucleus is equal to the average electron 
density in the plasma 0. In a recent review of solar fusion cross sections ||, the absence 
of an analytical study of nonlinear screening effects was underlined, while at the same time 
it was suggested that electron degeneracy and non-linearities of the Debye-screening could 
produce corrections to Salperer's formula roughly of the order of a few percent for solar 
reaction. On the other hand, a numerical integration of the Poisson-Boltzmann equation for 
a mixture of electrons and ions ||, yielded non-linear corrections of the order of 1%. 

The aim of the present paper is to study analytically non-linear (higher order) corrections 
for stellar nuclear reactions and determine their effects on the solar neutrino fluxes. The 
paper is organized as follows: In unit II the fundamentals of screened thermonuclear reaction 
are briefly reviewed. The non-linear screening formalism is derived in unit III and the effects 
of non-linearities on the most probable energy of interaction and the cross section factor 
are clarified. In unit IV the derived formula are implemented for various stellar reactions 
while in unit V the sensitivity of the solar neutrino fluxes to such non-linear corrections 
is investigated. Finally the conclusions of this paper are given in unit VI. An appendix 
at the end elucidates some misconceptions about the effects of electron screening on the 
astrophysical factor and the most probable energy of interaction. 
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II. FUNDAMENTALS OF ELECTRON SCREENING IN THERMONUCLEAR 

REACTIONS 



In the stellar interior the nuclear fusion cross section a (E) can be written as a product of 
a penetration factor P (E) times a nuclear factor a nuc (E) . The penetration factor is actually 
the probability that two positively charged particles will tunnel through the Coulomb barrier 
that separates them: 

a(E) = a nuc (E)P(E) (1) 
For two bare nuclei of mass numbers Ai,and A 2 respectively, that barrier is given by 

E. - (2) 

where the two nuclei can be considered sharp-edged spheres so that 

R=1A(A\' 3 + J£' 3 ) fm (3) 

Ai, A 2 being the mass numbers of the fusing nuclei. 

If higher relative angular momenta are considered in the collisions then the height of the 
barrier is increased by 

Ei = "' + f 2 (4) 
2/ir 2 V ; 

where \x is the reduced mass . 

In the stellar plasma the kinetic energy of the interacting particles is determined by a 
Maxwell-Botzmann distribution of velocities corresponding to a thermal energy 

kT = 0.086T 6 keV (5) 

where T 6 the temperature in millions degrees Kelvin. In a semiclassical approach it can be 
shown that the two nuclei can only interact if 



l<2-10- d R ifm)y jAT 6 (6) 

where A is the reduced mass number. Note that all energies are assumed to be center of 
mass energies unless specified otherwise. By means of the above condition it is easy to show 
that in stellar fusion reactions between light nuclei s— wave interactions dominate. In such 
a case the interaction is adequately described by a single particle potential of the form: 

V(r) = V N (r) + V c (r) (7) 

where Vn (r) is the nuclear potential and V c (r) is a potential which describes the Coulomb 
interaction and is not necessarily a pure Coulomb potential. Gamow first showed , in 
connection with the problem of a— decay that for two bare nuclei of charge Z 1 and Z 2 moving 
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with relative velocity u the probability to penetrate the Coulomb barrier is proportional to 
the factor exp (— 27m) where n is the Sommerfeld parameter: 



n(E) = Z x Z 2 e 2 



2h 2 



The cross section of the nuclear fusion reaction is then given by: 

a, (E) = ^-T t (E) 



E 



where 



T t (E) = exp 



h 



R 



V c (r) + 



1(1 



- Edr 



and r c is the classical turning point given by 

1(1 + 1) 



V r c ) + 



E 



(8) 



(9) 



(10) 



(11) 



The lower limit R of the integral is the radius of the nuclear forces given by (|3D (assumed 
to be practically zero) 

In the framework of the nearly perfect ionized gas, the presence of the electron cloud 
around the nuclei increases the reaction rates over their laboratory analogs. This screening 
effects has been studied by many authors and the most popular potential is the Debye-Huckel 
(D-H) screened Coulomb potential given by 



VAr) 



a 



(12) 



where u = r/ro and U e = ZiZ2e 2 r D 1 is the screening potential energy i.e. the Coulomb 
energy of the two colliding atoms at a distance equal to the D-H radius given by: 



D 



Aire 2 



J2 z l n i + n e0?j 



(13) 



where rij is the number density of ions with charge Zi , n e is the average electron density 
and 9 e the electron degeneracy factor. 

At this stage we need to define the limits of the Debye-Huckel model which are actually 
the limits of the results of this paper. This model assumes a nearly perfect gas at low density 
where the average Coulomb energy (E c ) between two adjacent nuclei is much smaller than 
the thermal kinetic energy of the plasma 



(E c ) «kT 



(14) 



Therefore, in the calculations that follow we have to bear in mind that the derived non-linear 
corrections cannot be applied beyond the realm of validity of this condition. 
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III. NONLINEAR SCREENING FORMALISM 



Let us assume a screening potential of the form (|T2"D . The cross section for thermonuclear 
reactions can be written as : 



E 



where x — — . 

The quantity / (x) is 



Kx) 



1 1 



e x(l-u) _ 



o V u 



where x = x (E) is the solution of the equation: 



xe 



E 



and / (0) is : 



1(0) 



ldu 



o V u 



7i 

2 



By introducing a multiplicative corrective term £ (x) such that 



7T 



I(x) = I(p)t{x) = -t{x) 



we obtain 



2 3 



73x 6 



£ (x) = e~ x H 1 1 1 

s w V 2 16 48 1024 5120 49152 

The screened reaction rate between nuclei i and j is: 

r%=(l + 6 ij )- 1 n i n j (av) ae 



+ 



(15) 



(16) 



(17) 



(19) 



(20) 



(21) 



where rii,nj are the number densities of nuclei % and j respectively, ^ is the Kronecker 
delta, and the thermalized cross section per pair of particles is 




(kT)' 



3 r°° 



S (E) exp 



E 



An(E)I(x) 



dE 



(22) 



Note that in the present paper the superscripts: nos, wes, sc, ss, indicate respectively no- 
screening, weak-screening, non-linear screening and strong screening regimes. 
Therefore 



(avY 




(kT)-*f (E° c ) S(E)exp 



dE 



(23) 
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where the screening correction factor f which is now evaluated at the most probable energy 
of the screened interaction is given by: 



2 1 

00 cc 



35x 4 23x 5 U9x 6 



ln/o = Trn (E° c ) \x - — - — + — — + — — + O x 7 ) (24) 

J0 K ; V 8 12 512 640 24576 L \J v 1 

Following the method of the steepest decent ||, the most effective energy of interaction E^ c 
between the two screened nuclei is: 

+ =0 (25) 

which yields (see Appendix I): 



E=E% 



E s c =[2I(x)Z 1 Z 2 e 2 ] J^kT ] j (26) 

where x = x (Eq c ) is the solution of equation (|17D . 
Hence 

E° c = £ 2 / 3 ( x ) E™ s (27) 

On the other hand by using the formula for the most effective energy of interaction for 
two bare nuclei we obtain: 

E S Q C = 1.220 • (zlZ 2 2 AT%) 1/3 e 2/3 (x) keV (28) 

Eq. (|TTD now reads: 

&o (ATi) r D(fm) 
In the weak screening (wes) approximation, to first order in x, we obtain: 

72 rr2 Arp2\ X / 3 /i X 



E™ es ~ 1.220 ■ [Z(Z^ATi) ( 1 — — ) (30) 



which corrects the result of ref |T(J according to which the energy should be shifted by 
(1 — x) , instead. 

Moreover, in the weak screening approximation, the integral / (x) can be written: 

IWES {x) ~ \ v 1 " I) (31) 



and therefore one obtains Salpeter's weak screening correction [10]: 

jwes _ e nnx ^ 
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Note that the assumption made in ref [M about the classical turning point being equal for 



both the screened and unscreened cases is not necessary JO]. In fact the actual approxima- 
tion is 



tT s l-x) 



(33) 



Regarding the rate of thermonuclear reactions , it has been shown that in the stellar interior 
where the temperature is T degrees Kelvin and the density is p in gj cm 3 the nuclear reaction 
rate is:. 



where 



2.62 • 10 29 2 XiXj 2 , _, 

— — — : — - — Too e ffT e cm sec 



3E™ S 
kT 



(34) 



(35) 



and S e ff is expressed in keV-barns. The quantities X it j, Zij, Aij, are the mass fraction, 
charge and mass number respectively of the nucleus i,j . 
To first order in r _1 |12J : 



S.„ (E ) ~ S (0) 



E=0 



e + 3 ^t) + ... 



(36) 



Although in Eq.(|3J) the screening enhancement factor f is now evaluated at Eq c the 
effective astrophysical factor S e ff is always evaluated at Eq os as the screening effects have 
already been worked out of the integral so the correction assumed in ref. [10] is not applicable 
(see Appendix II) 



IV. RESULTS FOR VARIOUS STELLAR FUSION REACTIONS 

In the solar region of maximum energy production (R/R & = 0.09) , the temperature and 
the density are respectively T 6 = 13.5 and p = 93.3-^3. Using recently calculated isotopic 



abundances of the solar interior |L3| the average internuclear distance given by 



a = = SlOOOp- 1 ^™ (37) 

(47rpA^ ) 

is found to be a ~ 11244 fm and the Debye-Huckel radius is = 25719 fm, where the 
weak electron degeneracy has been taken into account (9 e ~ 0.92) . Moreover, the thermal 
kinetic energy is kT = 1.161 keV and is always higher than the average Coulomb energy 
between two adjacent ions according to condition (114]) . Therefore the potential of Eq.([l2|) 
can be used in order to study the screening corrections to solar fusion reactions. The results 
are depicted in Table 1 that follows: 
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Table I. The electron screening corrections for various solar fusion reactions. 



Fusion TccictioTi 


fwes 

Jo 


fSC 

J 




A f n 


E% os {keV) E s c (keV) E c (keV) U e (ke 


H l {p,e + u e )H 2 


1.049 


1.048 


0.030 


0.5 


5.624 


5.471 


514 


0.056 


He 3 ( 2 #e 4 ,7) Be 7 
He 3 ( 2 He 3 ,2p) He 4 


1.212 


1.210 


0.129 


1.7 


20.860 


20.786 


1358 


0.223 


1.212 


1.210 


0.135 




19.952 


19.877 


1426 


0.224 


Be 7 (p, 7) -B 8 


1.212 


1.210 


0.161 


1.5 


16.671 


16.596 


1412 


0.224 


iV 14 (p, 7) O 15 


1.400 


1.396 


0.332 


0.8 


24.735 


24.606 


2111 


0.391 


C 12 (p >7 ) N 13 


1.335 


1.331 


0.272 




22.238 


22.127 


1876 


0.348 



Table I. The electron screening corrections for various solar fusion reactions with the follow- 
ing notation : /^ es , /Q C ,the screening factors for the weak and non-linear screening respectively, 
Af a (%) the non- linear correction of the weak screening in percent, Eq OS , Eq C the most probable 
energy of interaction for the two cases, E c the coulomb barrier of the reaction and U e coulomb 
energy of two ions at distance Td- Moreover the available corrections of the numerical approach [|| 
are also shown as A/ n (%) . 

Note that the statement || that non-linearities of the Debye screening might produce a 
correction to Salpeter's formula of the order of a few percent can be now upgraded. According 
to the above results, as far as the Debye-Hiickel mean field potential is concerned, the actual 
deviation is less than 0.5% for all solar fusion reactions. 

A very interesting fact is that such non-linear corrections are largely due to the presence of 
Eq . (p9|) which supersedes its linear oversimplified counterpart of the weak screening regime: 

x = (38) 

rpnos \ ' 

E 

That means that for most stellar conditions instead using Eq.(|24]) we can safely neglect 
higher order terms in order to use: 

In f Q = ixn (E° c ) x (E S C ) (39) 



where x (Eq c ) is always the solution of Eq. 

Moreover, the usual assumption that the screening enhancement factors are indepen- 
dent of the isotope is not necessarily valid when nonlinear corrections are considered. This 
is also clear from Eq.(p9|) where x depends on the reduced mass number. However, as can 
be readily seen from table I, for typical solar conditions, isotopic dependence is indeed prac- 
tically negligible. On the other hand the Gamow peak energy, which is usually considered 
screening independent, in the non-linear regime shows a maximum variation of the order of 
0.5% . 

For purposes of illustration nonlinear screening corrections have also been calculated for 
the cases considered by Salpeter 0. Following the notation of table I we have obtained the 
following results: 

a) In a red dwarf (main sequence star cooler then the sun), with typical central conditions 
p = 100, Tq — 8, 9 e — 0.8, the screening corrections for the pp reaction are: 



S 



fwes f SC A/(%) E™ s (keV) E s Q c (keV) E c (keV) U e (keV) kT (keV) . . 

1.116 1.114 0.132 3.873 3.848 514 0.076 0.688 ^ ^ 

b) Salpeter also considers the reaction N 14 (p, 7) O 15 for the interior of the hotter main 
sequence stars for a combination: p = 122, T 6 = 11.6, 9 e = 0.85. By taking into account 
non-linear corrections the results are: 

jwes jsc a/(%) E™ s (keV) E° c (keV) E c (keV) U e (keV) kT (keV) 

1.648 1.636 0.685 22.356 22.193 2111 0.499 0.997 ^ ' 

For the same reaction in the center of Sirius where p = 80 , T 6 = 20, and 6 e ~ 1. the 
corrections are: 

f «,es f sc A/(%) E nos^ keV) E «> ( keV ) E c (keV) U e (keV) kT(keV) 

1.209 1.207 0.120 32.145 32.037 2111 0.327 1.720 1 } 

c) In the more luminous main sequence stars where Hydrogen has been exhausted, 
the first stage of the two-stage reaction 3ife 4 — ► C 12 + 7 is the temporary formation of 
Be 8 , through the reaction He 4 + He A — > B 8 . The screening corrections for such a collision 
between two alpha particles, for indicative conditions considered by Salpeter are: 

i) p = 10 4 , T 6 = 150, 6 e = 0.87 



fwes Af{%) E™ s (keV) E° c {keV) E c {keV) U e {keV) kT (keV) 

1.049 1.048 0.017 109.346 109.141 1295 0.619 12.9 



(43) 



Despite the high density , obviously condition flT4| ) still holds here , therefore Salpeter's 
weak screening approximation is still reliable, 
ii) p = 10 6 , T 6 = 150, 9 e ~ 

jwea fsc Af(%) E™ s (keV) E° c (keV) E c (keV) U e (keV) kT (keV) 

1.491 1.474 1.136 109.346 107.680 1295 5.540 12.9 ^ ' 

In the above case (ii) , the electron screening effect is strong enough for the weak screening 
approximation to become inaccurate and therefore nonlinear corrections become important. 
Unfortunately, even nonlinear corrections cannot redeem the inaccuracies of the Debye- 
Hiickel potential which is only applicable to fusion reactions where condition ( III] ) holds. On 
the other hand the strong screening formula [Bf 



f: s = exp 



/ \ i/3 
0.205 (f) T- 



[z x + z 2 f 3 - zl ri - zl /s 



(45) 



where p e is the mean molecular weight per electron, can be cautiously used here as screening 
is strong, though not strong enough to assume absolute validity . By applying formula ( flop 
we obtain /q s = 1.498 which is reasonably close to its linear and non-linear counterparts. 
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V. EFFECTS ON SOLAR NEUTRINO FLUXES 



It has been suggested [14| that the electron screening effect may be important to the solar 



neutrino production as it enhances the reaction rates in the interior of the sun. More pre- 
cisely, if we consider the Gamow peak screening-independent then any variations of screening 
reflects on the cross section S (E) which is in fact multiplied by this screening factor so that: 

S sc (E) = f S nos (E) (46) 

In fact, it has been shown that the solar neutrino fluxes $ can be given as a function of the 
screening factors of the pp and CNO chains 0]. In that work the screening factor / was 
assumed to be isotope independent. Although in typical solar conditions this assumption 
was shown to be valid (see table I) when nonlinear screening effects are considered the 
derived formula should be modified as follows: 
Be 7 (e~,u e )Li 7 : 



Be 7 (p, 7 ) B 8 (e+ u e ) B 



8* 



,1V"'-"// \ -23.6/8 (f s fHe 3 +He i , K , 

= ®B Up+p) Up+ 7 Be) 7= (48J 

jHe 3 +He 3 



On the other hand for the CNO cycle which is governed by the slowest reaction 
A 14 (p, 7) O 15 the enhancement on the neutrino fluxes are 0: 
A 13 (e + v e ) C 13 and O 15 (e+, u e ) A 15 : 

3£,o = *3£> (f P+P r 22/S (U + ^n) (49) 
Conservation of luminosity yields the pp neutrino fluxes by means of 

%l + &Be + ®N + ®0 = %°P + ®B°: + ®T + $ 0° S (50) 

whereas the pep neutrino fluxes can be obtained by the observation that in any solar model: 

^pep _ &pp 

(hnos (jjnos ^ ' 

pep PP 

Instead of solving algebraic equations, an easier way of obtaining corrections to pp and pep 
neutrino fluxes is by using the proportionality formula [jT5[|: 

^ q0.14q0.03 q-0.06 /r \ 

^pp ~ iJpp o He 3 +He 3 He 3 +He 4 yoz) 

which readily yields: 

^pp _ j0.14j-0.03 f-0.06 /eo\ 

— Jp+pJHe : i+He3JHe3+He 4 \ 00 ) 

PP 

By using the above formulae the corrections to the solar neutrino fluxes both in the weak 
screening and non-linear screening regimes are shown in table II. 



10 



Table II. The nonlinear screening corrections of solar neutrino fluxes 



Neutrino source 


Qwes 1 (^nos 




H 1 (p, e + v e ) H 2 


1.000 


1.000 


H 1 [pe-,v e )H 2 


1.001 


1.001 


Be 7 (e-,u e )Li 7 


1.037 


1.037 


B 8 (e+,u e )B 8 * 


1.158 


1.159 


N 13 (e+z/) C 13 


1.227 


1.226 


O 15 (e+,z/ e ) A 15 


1.227 


1.226 



Admittedly, it has been argued [HJ that the weak-screening factors are more appropriate 
for the pp reaction rate where condition (|14D is fully satisfied, whereas for reactions with 
Z1Z2 > 4 Mitler's formula j|] should be used, instead. In our approach (Table II) we have 
used the weak-screening formula for all solar reactions following the recent suggestion of 
ref. Note that in ref. |2j, where Mitler's screening factors were used, the ratio $ sc /$ nos 
for the neutrino was found to be 0.995 while for the A 13 neutrino it was 1.13. This 
discrepancy is unimportant to the results of this paper according to which there is a negligible 
contribution of the investigated non-linear screening effects to solar neutrino fluxes well below 
the experimental errors of any currently imaginable neutrino detector. 



VI. CONCLUSIONS 

In this paper we have studied the non-linear effects of electron screening on stellar nuclear 
fusion rates calculating the respective corrections analytically. The formalism employed has 
been based on the Debye-Hiickel model. In typical solar conditions such non-linear effects 
are shown to be negligible proving Salpeter's linear approach to be sufficient for the study 
of solar nuclear reactions. Regarding the solar neutrino problem it was also shown that, 
non-linear screening leads to a negligible correction of the solar neutrino fluxes of the pp and 
CNO chains. Moreover, non-linear corrections are shown to be of some importance only in 
the intermediate screening regime, where the average Coulomb energy begins to challenge 
the average thermal kinetic energy and the Debye-Hiickel model begins to break down. 
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APPENDIX 

I. The most effective energy of interaction. 

The method of the steepest decent yields the maximum of the integrand of Eq. (f2~2|) : 

+ An (E) I (x)) = (54) 
dE \kT J e=e s c 

In fact the quantity x is energy dependent x = x (E) by means of Eq..(|l7]), therefore upon 
differentiation we obtain: 
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kT 



+ 2ir 



'dn(EY 
dE 



i «) + n (E° 



E=E, 



'dx (EY 
dE 







E=E, 



dx 







(55) 



where n (E) is given by (BI) and x s c is the solution of Eq..([D]) . If we assume that throughout 
the integral which appears in the thermalized cross section in Eq.(|22|) there is only a negligible 
variation of £ (x) then the above integral yields Eq.flSEp .This assumption is better than 
assuming that £ (x) = 1, which practically yields Eq os = E^ c , as the latter disregards a 
priori all the corrections of the screened interaction to the effective energy of interaction . 
II. Screening independence of S e ff . 

The basic formula for the thermalized cross section of the non-resonant screened ther- 
monuclear reaction is : 



(av) s 




S (E) exp 



~-4n{Ef)I(x) 



dE 



(56) 



We have to work out the screening factor f first otherwise the introduction of an S^t for 
the screened case replaces the integral above by an average (corrected) expression as follows 



12 



S (E) exp 



dE 



2E S C ( - 



1/2 



(57) 



Once the integral has been replaced the screening effect cannot be parametrized in the usual 
way by means of the screening factor f . 

Therefore, after we work out the screening enhancement factor f we obtain: 



(avY 




fo (E S C 



S (E) exp 



dE 



(5f 



where it should be noted that the evaluation of fo is performed at the most probable energy 
of the screened interaction. The remaining integral is actually unaware of the screening 
effects as the former is calculated by the method of the steepest decent around its maximum 
value, which is the (unscreened) most effective energy of interaction Eq os . 
Therefore, the energy E appearing in the formula: 



S eff = S(E ){l + r~ 



5 

12 



5 S' (E )E 
2 S(E ) 



S" (E ) E { 



S(E ) 



+ 



_— 2 



(59) 



1 /3 

can only be the quantity Eq os = 1.220- (ZfZ^ATg) , which corresponds to the no-screening 
regime. As a result, no screening correction can be incorporated into the frequently used 
formula: 



S eff ~S(0)\l 



5 S' (0) 
127 + S{0) 



E + -kT 
U 36 



4- S " (0) F 




(60) 



because Eq cannot be replaced by Eq c . Note that even if that was the case one should have 



to modify r = r (Eq os ) as well. 
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